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ABSTRACT 


The  nonlinear  development  of  the  negative-mass  instability 
is  investigated  for  a thin  E-layer  (A<<Rq)  located  inside  a 
cylindrical  waveguide.  A quasilinear  analysis  shows  that 
the  instability  saturates  by  the  development  of  an  energy 
spread  and  a corresponding  increase  in  thickness  A of  the  E- 
layer.  The  efficiency  of  microwave  generation  by  this 
instability  is  also  calculated. 
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1 . INTRODUCTION  AND  ASSUMPTIONS 


2 
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The  negative-mass  instability  has  received  considerable 

2-5 

attention  as  a mechanism  for  intense  microwave  generation. 

For  the  most  part,  however,  analyses  of  this  instability  have  been 
limited  to  the  linear  regime.  The  purpose  of  this  paper  is  to 
investigate  the  nonlinear  development  of  the  negative-mass  instability, 
with  particular  emphasis  on  the  implications  for  intense  microwave 
generation.  The  analysis  is  carried  out  for  a thin  (A«Rq)  , 
relativistic  nonneutral  E-layer  aligned  parallel  to  a uniform  axial 
magnetic  field  and  located  inside  a grounded  cylindrical 

waveguide.  The  nonlinear  consequences  of  the  negative-mass  instability 
are  investigated  within  the  framework  of  a simple  quasilinear  model 
that  makes  use  of  total  energy  balance  for  the  system.  One  of  the 
most  important  features  of  the  analysis  is  that  the  instability 
saturates  by  the  development  of  an  energy  spread  and  a corresponding 
increase  in  thickness  A of  the  E-layer.  The  corresponding  efficiency 
of  microwave  generation  is  also  calculated. 

As  illustrated  in  Fig.  1,  the  present  analysis  is  carried  out  for 
an  unneutralized  E-layer  aligned  parallel  to  a uniform  external  magnetic 

A 

field  Bq£z  and  located  inside  a cylindrical  conducting  waveguide 
with  radius  Rc.  Cylindrical  polar  coordinates  (r,6,z)  are 
Introduced,  and  the  main  assumptions  are  enumerated  below. 

(a)  The  E-layer  is  thin,  i.e.,  A/Rq«1,  where  A*R2~R^  is  the 
radial  thickness  and  Rq“(R1+R2>/2  is  the  mean  radius  of  the  E-layer. 

(b)  The  electron  motion  is  ultrarelativistic  (yq»1)  and  the 
mean  equilibrium  motion  of  the  E-layer  is  in  the  azimuthal 
direction  (i.e.,  V^-0,  where  is  the  mean  axial  velocity  of  an 


electron  fluid  element). 
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(c)  It  Is  further  assumed  that  v/yq  <<;1>  where  v*Nge  /me 

is  Budker's  parameter,  Ng  is  the  number  of  electrons  per  unit  axial 

2 

length  of  the  E-layer,  c is  the  speed  of  light  in  vacuo,  Ygmc  i® 
the  azimuthal  electron  energy  at  equilibrium  radius  Rq,  and  -e  and 
m are  the  charge  and  rest  mass,  respectively,  of  the  electron. 


2 . THEORETICAL  MODEL 


The  equilibrium  and  negative-mass  stability  properties  have  been 

investigated  in  Refs.  7-9  for  the  choice  of  a sharp-boundary  equilibrium 

in  which  the  electrons  have  a rectangular  density  profile,  i.e.,  n^(r)= 

Ne/(27rRQA)=const. , for  R^r^,  and  n®(r)=0,  otherwise.  Here  R2 

is  the  outer  radius  of  the  E-layer  and  R^=R2~A  is  the  inner  radius. 

For  a thin  E-layer,  the  resulting  dispersion  relation  for  the  negative- 

7-9 

mass  instability  can  be  expressed  as 


where  w is  the  complex  eigenfrequency , u)c=eBg/YQmc  is  the  electron 
cyclotron  frequency,  H is  the  azimuthal  harmonic  number,  and  use  has 

8 

been  made  of  the  cold-fluid  limit  of  the  kinetic  dispersion  relation 
for  the  negative-mass  instability.  In  Eq.  (1),  the  geometric  factor 
g is  defined  by  g=J,^/(b_+b+)-k^RQ/(d_+d+)  , and 


b.*=  - 


UJl(pR2)N’(pRc)-j;(pRc)N£(pR2) 


+ R2pJ^(pR2)N^(pRc)-J^(pRc)N'(pR2)  ’ 


R2PJi(pR2)N£(pRc)_JJl(pRc)N£(pR2) 
V ‘ Ujl(pR2)NJl(pRc)-JJl(pRc)NJl(pR2)  * 


b_=UJl(pR1)  / [RjpJ'  (pRx)  ] , d_«l/b_. 


where  J^(x)  and  N^(x)  are  Bessel  functions  of  the  first  and  second 
kind,  respectively,  the  prime  (')  denotes  (l/p)(d/dr),  and  p is  defined 


<* 


222  1/2 

by  p*(w  /c  -k  ) , where  k is  the  axial  wavenumber. 

The  nonlinear  consequences  of  the  negative-mass  instability 

can  be  investigated  within  the  framework  of  the  conservation 

equation  for  net  energy  balance.^  We  Fourier  decompose 

perturbed  quantities  according  to  6A(£,t)*=  £ 6A^  k(r,t)exp(i£0+ikz) , 

H,k  ’ 

where  k=  2irq/L  and  q is  an  integer.  Energy  conservation  can  then  be 
expressed  as 

..  fR  ^ 


It  £2u  fc  drre£,k(r't)  +h  K=0- 

k 0 


where  z ^ k*|6JgA  k(r,t)  1 2/8n+|6^  k(r,t)|2/8ir  is  the  perturbed  field 
energy  density,  and 


fL/2  r 2i 

T rR^  r 

dz  < 

16  drr 

o 

CM 

X 

1 

— * 

is  the  average  electron  kinetic  energy  per  unit  length.  Here  L is  an 
effective  periodicity  length.  Within  the  context  of  a simple  quasilinear 
model,  the  spectral  energy  density  evolves  according  to^ 


where  the  growth  rate  u>^’  =Imoj  is  determined  adiabatically  in  terms 
of  other  system  parameters  from  the  linear  dispersion  relation  in 
Eq.  (1).  Moreover,  for  a thin  E-layer  with  V^»0  and  near-laminar 
flow  in  the  azimuthal  direction,  we  approximate  the  electron  kinetic 
energy  per  unit  length  by^ 


t R „ j 

K*2ir  Cdrrn(r,t)  (y-l)mc  , 
Jn  e 


where  n®(r,t)=/^^  ^7  ne^»£^  is  the  avera8e  electron  density, 

2 2* 2 2 ] /2  2 

and  ymc  = (l+r  u>*/c  ) x/ me  is  the  azimuthal  energy  associated  with  the 

Q A a 

average  electron  motion  (V0=ruic/y).  Here  ii>c“eBQ/mc  is  the  nonrelativlstic 
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electron  cyclotron  frequency. 

Substituting  Eqs.  (4)  and  (5)  into  Eq.  (3),  and  assuming  that 
the  electron  density  profile  maintains  a rectangular  shape  with 
variable  thickness  A(t),  we  obtain 

Nemc2F(A)  - l <4,kVk(0>exP{2j  dt'a)J’k(t')  } , (6) 

! , k 0 


where 


rR 

W , (0)H2ir  cdrre#  . (r,t*0)  characterizes  the  energy  of  the  initial 
X ,1c  J Q X , K 


field  perturbations,  Ng  is  the  number  of  electrons  per  unit  axial 

length,  and  F(A)>0  is  defined  by 

2,3/2  2,3/2.  /,^l/2\‘ 


( 2x  [(l+x?)3/2-(l+x2)3/2)  3x,(l+x2>1/2  i i 

1 - 2 272  f ’ <?) 

\ X2  x^  / 


with  x0  = ui  R„/c  and  x.-w  (R„-A)/c.  Consistent  with  experimental 
2 c 2 1 C 2 

observations,  in  obtaining  the  left-hand  side  of  Eq.  (6)  we  have 
approximated  dR^/dt^O  but  retained  terms  proportional  to  dA/dt. 


3.  NONLINEAR  EVOLUTION 


For  g>0,  it  is  evident  from  Eq.  (1)  that  the  linear  growth  rate 

i k 

o)1*  vanishes  once  the  layer  thickness  has  increased  to  the  critical 


value 

h-T-  <¥*>m  ■ <«> 

f lit)  lyn 
c '0 

Therefore,  making  use  of  Eq.  (6),  the  quasilinear  evolution  of  the 
system  can  be  summarized  as  follows.  At  t»0,  the  initial  thickness 
of  the  E-layer  is  A^»A(t*0).  As  the  field  perturbations  grow,  A(t) 
continues  to  increase  fdA/dt>0  since  F(A)>0]  until  A(t-x«>)-Aj  and  the 
instability  ceases.  (Here  we  assume  that  a single  l value  is  excited.) 
Concomitant  with  the  increase  in  E-layer  thickness  is  an  increase  in 


azimuthal  energy  spread  of  the  electrons  composing  the  layer. 


6 


I 

! 

i 

! 


The  efficiency  n of  microwave  generation  is  defined  as  the  ratio 

of  electromagnetic  field  energy  generated  by  the  instability  to  the 

initial  kinetic  energy  of  the  E-layer.  In  the  ultrarelativistic 
2 2 

limit  with  Yq>>1,  X2»1,x^»1  and  A/Rq«1,  some  straightforward  algebra 
that  makes  use  of  Eqs.  (6)  and  (7)  shows  that  o can  be  approximated  by 


n = 2R^  (W 


(9) 


For  future  reference,  we  also  introduce  the  average  power  per  unit 
axial  length  associated  with  the  microwave  production,  i.e., 

?*.k  -T  (10) 

where  T is  the  duration  of  the  radiation. 

Equations  (1)  and  (6)  have  been  solved  numerically  for  the  self- 

consistent  nonlinear  evolution  of  the  E-layer  thickness  A(t).  For 

present  purposes,  we  consider  transverse  electric  (TE  mode)  perturbations 
2 2 2 2 2 2 

with  k Rn=A  -a.  (R„/R  ) and  J'(a„  )=0.  Typical  results  are  summarized 
0 An  0 c A An 

in  Figs.  2 and  3 for  Rc/Rq=1.5  and  Yq=4.5.  In  Fig.  2,  we  plot  normalized 

E-layer  thickness  A(t)/Rg  versus  wct  for  several  different  values  of 

Budker's  parameter  v=Nge  /me  . Evidently,  the  time-asymptotic 

thickness  of  the  E-layer  (A^)  is  a slowly  increasing  function  of  v 

[see  also  Eq.  (8)].  On  the  other  hand,  the  time  scale  required  for 

saturation  decreases  rapidly  with  increasing  v,  which  is  a consequence 

£ k 

of  the  fact  that  the  growth  rate  increases  with  v during  the 

initial  stages  of  instability  when  A«Af  [Eq.  (1)]. 

We  note  from  Eqs.  (8)  and  (9)  that  the  efficiency  n of  microwave 


| 


generation  is  also  an  increasing  function  of  Budker's  parameter  v. 

To  illustrate  this  point,  the  normalized  average  microwave  power  per 
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unit  length  P/u)c(mc  /e)  is  plotted  versus  v in  Fig.  3 for  Rc/Rq=1.5, 
Yq=4.5,  and  several  different  waveguide  inodes  (£,n).  As  an  example, 
for  the  mode>  t*ie  average  microwave  power  per  unit  length 

P increases  five-fold  when  v is  increased  from  0.07  to  0.15.  Also 
shown  in  Fig.  3 is  a plot  of  efficiency  n versus  v. 

4.  CONCLUSIONS 

Several  important  conclusions  follow  from  the  present  analysis. 

For  example,  the  time  evolution  of  the  E— layer  thickness  exhibits  a 

very  sensitive  dependence  on  Budker's  parameter  v [Fig.  2].  Moreover, 

the  average  power  and  efficiency  of  microwave  generation  are  enhanced 

considerably  by  an  increase  in  electron  density  [Fig.  3].  Introducing 

an  axial  energy  spread  into  the  present  analysis  reduces  the  growth 
0 k 59 

rate  w*’  substantially,  ’ thereby  slowing  the  rate  at  which  the 
E-layer  thickness  increases.  This  behavior  is  consistent  with  recent 
experimental  observations  by  Destler  et  al.  (cf . , Figs.  8 and  11  in 
Ref.  3). 
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FIGURE  CAPTIONS 

Fig.  1 Equilibrium  configuration  and  electron  density  profile. 

Fig.  2 Plot  of  A(t) /Rq  versus  a>ct  for  Rc/Rq=1.5,  Yq=4.5,  Aj/^0= 

-2  - 2-3 

2x10  , ^(0)/(Nemc  )=8xl0  , and  several  different  values 

of  Budker's  parameter  v. 

2 2 

Fig.  3 Plot  of  normalized  microwave  power  per  unit  length  P/wc(mc  /e) 

versus  v [Eqs.  (10)  and  (6)]  for  Rc/Rq=1.5,  Yq=4.5,  A^/Rq=2x10 
- -3  2 

W . (0)=8xl0  N me  , and  several  different  modes  (£,n). 

JC  K.  0 

Also  shown  is  a plot  of  efficiency  n versus  v. 


